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ABSTRACT 

We study the evolution of the stellar and dark matter components in a galaxy cluster 
of lO'^^M© from z = 3 to the present epoch using the high-resolution collisionless 
simulations of Ruszkowski & Springel (2009). At z = 3 the dominant progenitor halos 
were populated with spherical model galaxies with and without accounting for adia- 
batic contraction. We apply a weighting scheme which allows us to change the relative 
amount of dark and stellar material assigned to each simulation particle in order to 
produce luminous properties which agree better with abundance matching arguments 
and observed bulge sizes at z = 3. This permits the study of the effect of initial com- 
pactness on the evolution of the mass-size relation. We find that for more compact 
initial stellar distributions the size of the final Brightest Cluster Galaxy grows with 
mass according to r cx Af^, whereas for more extended initial distributions, r cx M. 
Our results show that collisionless mergers in a cosmological context can reduce the 
strength of inner dark matter cusps with changes in logarithmic slope of 0.3 to 0.5 
at fixed radius. Shallow cusps such as those found recently in several strong lensing 
clusters thus do not necessarily conflict with CDM, but may rather reflect on the ini- 
tial structure of the progenitor galaxies, which was shaped at high redshift by their 
formation process. 
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1 INTRODUCTION 

Brightest cluster galaxies (BCGs) are the most massive 
elliptical galaxies in the Universe at the extreme end of 
the galaxy luminosity fu nction and perhaps "special" (see 
iParaniape fc Shetbll20lih . These objects are intriguing be- 
cause they do not follow the same scaling relations as nor- 
mal gia nt ellipticals. Using the Slo an Digital Sky Survey 
(SDSSI. Ironder Linden et al.l ll2007f) found that BC Gs devi- 
ate systematically from the lFaber fc JacksonI (| 19761 ) and the 
iKormendvl (| 19771 ) relations with lower velocity dispersions 
and larger sizes respectively. Early theoretical studies inves- 
tigated the role of cooling flows (Fabian 1994) in BCG for- 
mation. This hypothesis is disfavoured by observations with 
Chandra and XMM-Newton which show that such flows are 
much weaker than required, as are the star formation rates 
in the ce ntral galaxies of clusters (jMcN amara et aL I I2OOOI : 
iFabian et al. 2001). Another scenario has BCGs growing by 
feeding on smaller galaxies (minor mergers). This is the 



notion of galactic canniba lism (|Ostriker fc HausmanI [l977l : 
iHausman fc Ostrikeij[l978l ). 

In the ACDM scenario, structures build-up hierarchi- 
cally through accretion and mergers of smaller progenitors. 
Groups form before clusters and have sufficiently low rela- 
tive velocities that galaxy-galaxy mergers can occur before 
cluster formati on, thus enhanc ing the formation of a massive 
central galaxy. iDe Lucia fc Bl aizot ( 20071) use semi-analytic 
galaxy formation models to show that in the ACDM cosmol- 
ogy, BCGs form primarily through in-situ star-formation at 
high redshifts, 2^3, with subsequent mass growth dom- 
inated by non-dissipational merging. Similar processes are 
seen in hydrodynamic al cosmological simulations of mas- 
sive galaxy formation (|Naab et al.l l2009l : lOser et al] I2OIII : 
iFeldmann et al.ll201ll ) although these are typically less effec- 
tive at suppressing star formation at late times than semi- 
analytic models, leading to galaxies which are "younger" 
than those observed. 

From this, it seems that many aspects of the late assem- 
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bly of BCGs can be modelled without considering the early 
star formation phase. .Dubinskj (1998.) was an early exam- 
ple of such work based on numerical N-body simulations of 
collisionless mergers of galaxies in a cosmological context. 
Dubinski found that the properties of cluster BCGs can be 
naturally explained by merging of galaxies which have al- 
ready formed their s tars at high redshift. Th i s idea was fur- 
ther investigated by iRuszkowski fc Springel (|2009l ') (RS09) 
who studied the deviations of BCGs from the Kormendy 
and FJ relations in a ACDM simulation. 

In relaxed clusters, BCGs reside at the bottom of the 
potential well, making them ideal probes of the distribution 
of dark ma tter from kpc to Mpc scales. ISand et al.l (|2002l . 
|2004 [ 20081 ) studied a selection of clusters combining stel- 
lar dynamical modelling and strong gravitational lensing in 
order to infer the inner slope of the dark matter density pro- 
files. Their studies found values for the logarithmic slope of 
the dark matter density profile 7 = — ciln(p)/dlnr < 1, 
at odds with the predictions of dark-matter-only simula- 
tions of halo formatio n which generally f ollow the NFW 
profile with 7n f w = 1 (|Navarro et al.|[l997h . More recently, 
Newman et al. (|201lh revisited the study of Abell 383 by 



Sand et al.l ( 2008l l. combining stellar kinematics, strong and 



weak lensing and X-ray data to deduce an inner slope of 
7 = CSQlg gg at 95 percent confidence. These authors sug- 
gest this may indicate a genuine problem with our under- 
standing either of baryonic evolution or of the nature of the 
dark matter. 

Indeed, some recent hydrodynamical simulations of 
galaxy formation in clusters indicate a steepening of the 
inner - slope of the dark matter profile ( Gnedin et al.l |2004| . 



I2011I : [s ommer-Larsen fc LimousinI l20ld ). exacerbating the 
core-cusp proble m . More gentle co ntraction is seen in other 



simulations ( e.g. Duffv et al 



vestigated bv lNewman et al 



I2OIGI) but not on the scales in- 
' (|201ll). A simulation in which 



dynamical friction significantly erodes the cusp a t the c enter 
of ellipticals was presented by [Johansson et al.l (|2009l) but 
on the scale of the central galaxy of a small group. How- 
ever, the final masses of the central galaxies in these simu- 
lations are in general a factor of two or three h igher than 
expected from abundance matching argumen ts (|Guo et al.l 
l2010l : lMoster et al.|[2010l : lBehroozi et alfeoi O^. implying the 
need for a signicantly improv ed treatment of baryonic astro- 
physics. lEl-Zant et al.l (|2004l ') claimed that shallower cusps 
could be produced in a cluster through dynamical heating 
by the galaxies. However, they treated galaxies as unstrip- 
pable point masses which is too unrealistic to address the 
issue in quantitative detail. 

At this stage it still seems interest ing to address the 
second question bv lNewman et al.l l|201ll ): is the presence of 
shallow dark matter cusps at the centre of clusters a signifi- 
cant challenge to CDM? We use the RS09 simulation to test 
whether such cusps can be created through dry (i.e., gas- 
free) mergers. Recent observations of massive ellipticals at 
z = 2 have shown that they wer e more compact than sim - 
ilar mass galaxies today (see e.g. Ivan Dokkum et al.|[2008l ). 
Dry, predominantly minor mergers have also been proposed 
as a possible mechanism to drive the required size evolution 
(e.g., iNaab et al.ll2009l : iBezanson et~alll2009l ) 

In this context a significant limitation of the RS09 sim- 
ulations was that the galaxies they inserted at z = 3 had 
stellar masses an order of magnitude larger than expected 



from abundance matching arguments l|Moster et al.l |201G| ) 
and were assumed to follow the present-day mass-size rela- 
tion. Here we remedy the inconsistencies between the simu- 
lations and observations by using a method that re-assigns 
the mixture of stellar and dark matter in each simulation 
particle. This enables us to study the evolution of stellar 
and dark matter distributions for different levels of initial 
compactness and stellar mass. 

In Section 2, we give a description of the simulations 
as well as of the weighting scheme used in this study. We 
also present results on the mass and size growth of the BCG 
for different initial assignments of stars and dark matter. 
In Section 3, we look at how the initial slope of the dark 
matter evolves from z — 3 to the present. We discuss our 
results and conclude in Section 4. 



2 NUMERICAL METHODS 

2.1 Simulation 

" " ^dy 

tail in IRuszkowski fc Springeil (|2009l '). (RS09) and we give 
only a short summary here. A cluster mass dark matter 
halo of 10^^ Mo was identified in the Millennium Simulation 
(jSpringel et al.ll2005l ). The cosmological parameters of this 
simulation are Q.m = 0.25, Q.a = 0.75, a scale-invariant slope 
of the power spectrum of primordial fluctuations (n = 1.0), 
a fluctuation normalization og, = 0.9, and a Hubble constant 
Ho = lOO/ikms"^ Mpc~^ = 73kms-^Mpc-^ 

The cluster was then re-simulated using a zooming tech- 
nique with a mass resolution of m = 1.57 x 10^ Mq and 
comoving softening length e — 2.0 kpc. At redshift z = 
3, the 50 most massive progenitors of the flnal cluster were 
identified and replaced by spherical equil ibrium models in 
which stars were distributed according to a lHernauistl (Il990l) 
profile embedded in a NFW dark matter halo ( Navarro et al.l 
[l997h . The masses of the dark matter and star particles were 
set to be identical rudm = ni, = 1.57 x 10*^ Mq and they 
were assigned a softening length of e = lh~^ kpc. comoving, 
half that of the other dark matter particles from the original 
resimulation. Throughout this paper we use h — 0.73 and 
our mass and length units are thus in kpc and Mq. 

Two simulations were run with different initial galaxy 
models: one in which th e dark matter w as adiabati- 
cally contracted following iBlumenthal et al] (jl986) and 
one where it retained an undisturbed NFW profile. We 
shall refer to these a s mo dels A and B respectively. 
The IBlumenthal et al.1 (119861 ) formalism over-predicts the 
amount of contraction observed in many hydrodynamical 
simulations (Giiedin et ah 2004), however through the in- 
clusion of uncompressed and compressed dark halo models, 
we can probe two alternative regimes. If dry-merging is in- 
deed the main driver in the late assembly of BCGs and if the 
RS09 initial galaxies were realistic, then the "real" solution 
would lie between these two models. 

In fact, however, the galaxies which RS09 inserted at 
reds hift z = 3 followe d the present-day mass-size relation 
from lShen et~aLl (|2003h and assumed a stellar to dark mat- 
ter mass ratio m^^/M — 0.1. This value is too large by 
a factor of 10 according to recent results from matching 
the observed high-redshift abundance of massive galaxies 
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l|Moster et al.ll2010l : iBehroozi et"alll2010l : IWake et al.ll201ll ) . 
We note that many other simulations investigating similar 
processes hav e like-wise assumed over-massive stella r com- 
ponents fe.g- lNipoti et al.ll2009l : iRudick et aLlbOllI ). How- 
ever, we point out that for a given mass resolution, more 
massive stellar bulges are represented by a larger number of 
particles which considerably improves the numerical conver- 
gence of the simulations. 

Not surprisingly, the final merger remnants in RS09 
were also too massive m* ~ 10^"^ Mq with half-mass radii 
which were too large (~ 100 kpc) compared to real BCGs. 
Thes e generally do not exceed re ~ 50kpc (|Bernardi et alJ 
I2OO7I ). 

In order to study more consistently the change in the 
slope of dark matter density profiles at the centre of ACDM 
clusters, we need to address the question with galaxies that 
have stellar properties consistent with observations. We also 
need to test whether the final stellar mass of the merger rem- 
nant a grees with t he stellar to halo m ass relati on (SHM) at 
z = (|Guo et al.l :2010: Mo ster et al. 2010: Behroozi et all 
I2OIOI ). We employ a weighting procedure to re-assign the lu- 
minous component of every initial galaxy such that only one 
percent of its total dark matter mass is locked in stars. For 
the range of halo masses (10^"^ M© — 10^^ Mq) that we pop- 
ulate, abundance matching implies that this ratio is rather 
constant. Additionally, our weighting scheme enables us to 
change the sizes of the luminous components to study the as- 
sembly of the BCG for different levels of compactness while 
keeping the total initial stellar mass fixed. We present this 
scheme in the next section. 




E [kmVs^] 



2.2 Weighting Scheme 

Here we describe our method for re-assigning simulation 
particles in the initial conditions so that the light to stel- 
lar mass ratio and the size of the luminous component of 
a galaxy can be varied. O ur scheme is similar to that of 
iBullock fc JohnstonI (|2005l ) and is applicable to spherically 
symmetric density distributions. A galaxy is represented by 
a distribution of A'^ particles of mass Mtot/N = irip and 
phase-space coordinates (x, v) which generate the potential 

"1> = 'I'dark + ^stellar = "I'NFW + "I>Hcrnquist ■ Each particle of 

energy E = + $ now simultaneously represents dark 
matter and stars in different amounts according to a weight 
function iu{E) — ft,{E)/ f{E). This is the ratio between the 
stellar distribution function (df) and the total df. One can 
construct the stellar df in the following way. 

In computing /,(_£), one assigns the particles a spheri- 
cal number density distribution v{r) and solves the Edding- 
ton inversion formula: 



(fv{^) 1 dv 



(1) 



where * = -$ -|- $0 and £ = -E -I- $0 = 5' - /2 are 
the relative potential and total energies respectively. The 
potential used to generate the model galaxies in RS09 was 
a linear combination of a Hernquist and an NEW potential, 
both of which tend to zero in the limit r goes to infinity thus 
<l?o = 0. We choose v[r) to follow a Hernquist profile (1990): 



v{r) = 



r(r -f aY 



(2) 



Figure 1. Differential energy distribution for the proto-BCG for 
the total mass and light contributions for three test cases = 
^'iJRSog/-'-' '^211309/2, J'eRsog /3, rengog/5. Although the rej^gog/S 
histogram intersects the total differential energy distribution, the 
particles in those energy bins are below our spatial resolution of 
e = 1 kpc. 



The total df is the ratio of the differential energy dis- 
tribution N{E) = dM/dE = f(E)g{E) and the density of 
states g{E) which is solely defined by the potential $: 



g{E) = {47T'f / rV2(i5 - <E-(r))dr. 



(3) 



Since $(r) and N{E) can both be measured directly from 
the simulation's initial condition, this determines f{E). 

Note that the only free parameter we have introduced 
is the scale radius a for the light distribution which is re- 
lated to the half-light radius by a = re/(-\/2 -|- 1). Within 
certain limits we can vary the relative mixture of dark and 
luminous matter to represent, for example, less massive and 
more compact bulges. This has the advantage of allowing us 
to study aspects of the assembly of massive galaxies from 
2: = 3 to z = without having to run additional CPU in- 
tensive simulations, simply by tracking the weights to the 
final merger remnant. This scheme allows multiple interpre- 
tations of a single simulation. However, for our purposes the 
stellar mass within each galaxy is kept fixed, its value being 
dictated by abundance matching arguments, so it is the stel- 
lar and dark matter distributions which vary with the total 
mass distribution held fixed. Note that this implies that the 
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Figure 2. pr^ profiles vs. radius for the four different realisations. 
The target functional forms of each realisation are overplotted in 
solid colored lines. 

initial dark matter distributions in the adjusted galaxies are 
no longer those expected naturally in ACDM. Thus while 
we can address issues of how the mixing of the two compo- 
nents changes inner profile shapes and is affected by initial 
compactness, we cannot expect the final DM distributions 
to be realistic. 

Figure 1 shows the differential energy distribution for 
the most massive of the 2 = 3 galaxies (the dark halo and 
stars of the main BCG progenitor) as a function of energy. 
Overplotted are stellar differential energy distributions for 
four reinterpretations with different galaxy sizes: the orig- 
inal effective radius used in RS09 and the same reduced 
by factors of two, three and five. To show more explicitly 
that our method works we present pr^ profiles in Figure 2. 
Note that changing the sizes of the stellar component for 
a fixed m,/Mhaio ratio, implies substantial changes in the 
inner slope of the dark matter profile. The reduction of the 
stellar mass by a factor of 10 from that assumed in RS09 
also means that both the uncontracted and contracted mod- 
els now have overly concentrated dark matter distributions 
in the centres of the galaxy subhalos, except at very small 
radii where the reduced radii can lead to an increase of stel- 
lar density relative to RS09. 

The maximum extent to which we are able to rescale 
the stellar component is set by the total mass profile. This 
is saturated by the stars alone at the softening radius if the 
RS09 sizes are reduced by a factor of ~ 5 (see Figure 1). 

The initial and final light profile shapes in the RS09/1 



interpretation will be the same as in the original simulation, 
as only the stellar masses of every galaxy are changed (the 
inner dark matter profiles will differ, however, since they 
now contain the additional mass which used to be assigned 
to stars). This also m eans that the galaxies no longer lie on 
the lShen et"al] (120031 1 stellar mass-size relation. In order to 
put them back on it (within the scatter) we need to reduce 
the sizes by factors of ~ 5. 

Recent observations show, however, that z = 2 ellip- 
tical galaxies were more compact tha n implied by the lo- 
cal relation (|van Dokkum et al.l |2008| ) . Unfortunately, our 
spatial resolution limit does not permit us to consider such 
small sizes. We stress that these observations still need to 
be treated with caution as the galaxy stellar masses are es- 
timated photometrically. Martincz-Mariso et al. (2011) ar- 
gue that dynamical masses of compact galaxies at redshift 
z = 1 may be six times lower than some photometric esti- 
mates. Nevertheless if the photometrically determined stel- 
lar masses of galaxies at redshift z = 2 are even approx- 
imately correct, the galaxies should be even smaller than 
we assume in this paper. As we will see, the exercise pre- 
sented here can nonetheless give insight into the puffing-up 
of BCGs by minor mergers and its dependance on the initial 
compactness of the galaxies. 

2.3 Results for the BCG evolution 

2.3.1 Size growth of the BCG 

Fixing the stellar masses within all haloes according to abun- 
dance matching arguments (and hence reducing them by an 
order of magnitude from those originally assumed by RS09), 
we studied four assumptions for the compactness of the 
galaxies ~ i^brsoq /!> '"ensog/^i '"eRsog/3, each 
of our two simulations. The trends arevery similar in the 
two models. We find that the relative growth in size of the 
BCG from z = 3 to « = is not identical for different as- 
sumptions about initial concentration. This is illustrated in 
Table 1, where we give the growth factors as characterised 
by the increase in half-light radius. The size of the BCG at 
2; = in the original RS09 simulations was high by a factor of 
two compared to local BCGs. However, in the RS09/5 case, 
the sizes of the initial galaxies correspond to those of low- 
redshift largest ellipticals of similar stellar mass (within the 
scatter of the Shen et al. 2003 relation) and the most mas- 
sive BCG progenitor grows by a factor of 4, reaching a size 
that is in reasonable agreement with observed BCGs. This 
supports the idea that massive ellipticals can grow rapidly in 
size t hrough repeated minor mergers in a cosmological con- 
text (iNaab et al.ll2009l : lOser et al.li20ld. 1201 ll : IShankar etall 

Iml). 

2.3.2 Mass growth of the BCG 

Turning to the mass growth of the BCG, the main pro- 
genitor starts with a stellar mass of ~ 2.5 x 10^^ M0 at 
2 = 3 and the final merger remnant at 2; = has a to- 
tal stellar mass of 6 to 7 x IO^^Mq. For the cluster mass 
we are considering here (~ 10^^ M©), abunda nce match- 
ing suggests m,/Mhaio = 0.001 (see Fig. 2 of IGuo et al.l 
(2010)). Considering that semi-analytic models predict that 
80 percent of the stars which end up in a BCG are already 
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Run 


re/kpc 


ref /reiA 






Mf/A-hB 


RS09/1 


38 


2.7 


2.9 


2.5 


2.9 


RS09/2 


19 


3.2 


2.7 


3.0 


2.6 


RS09/3 


13 


3.6 


2.5 


3.1 


2.5 


RS09/5 


8 


4.3 


2.4 


3.8 


2.4 



Table 1. Initial half-light radii re^, size and mass growth factors 
(defined as the ratio of initial and final half-light radii and stellar 
masses) for the eight interpretations. The initial stellar mass of 
the largest BCG progenitor is 2.5 X IO^^Mq in all cases. 

formed hy z — 3 (|De Lucia fc BlaizotI |2007| ) our value of 
m»/M{z = 0) = 0.0006 is a considerable improvement over 
the value prior to rescaling (m,IM{z = 0) ~ 0.01). We 
note that this factor is still somewhat lower than expected 
from semi-analytic models. A quick query in the Millenium 
database indicates that only 30 percent of the stars in the 
z = BCG in the cluster come from the 50 most massive 
progenitor haloes at z = 3. The missing contribution comes 
from smaller subhaloes which were not populated with stars 
in the simulation. Figure 3 shows the normalised radial dis- 
tribution of simulation BCG stars coming from different 
galaxies at z = 0. The two panels illustrate the process of 
mass aggregation for two extreme interpretations (RS09/1 
and RS09/5). The stars coming from the proto-BCG (in- 
situ) are represented by the solid black line and those ac- 
creted from other galaxies by the dashed-dotted line. Con- 
tributions from individual galaxies are shown by coloured 
lines. Between z = ?, and z = the BCG was subjected 
to two mergers with mass ratios of about 3:1 and six with 
about 10:1. These are responsible for the inside out growth 
of the BCG. 

For the most compact galaxy interpretation (RS09/5), 
30 percent of the BCG mass comes from accretion which 
dominates beyond ~ 30kpc. This is in contrast with the 
extended galaxy (RS09/1) where this transition occurs at 
60kpc. Note that little of the accreted mass reaches the cen- 
tre. Most gets deposited on the outskirts leading to inside- 
out growth. 

The reason why the mass growth factor decreases with 
increasing compactness is that the initial compact galaxies 
have more stars on higher binding energy orbits than their 
less compact counterparts (as can be seen in Figure 1). These 
stars are harder to unbind. Thus, as encounters take place, 
less stellar mass is stripped and deposited on the BCG. More 
remains bound to the original galaxies. 

Another interesting aspect of our numerical experiment 
is the scaling between total stellar mass and effective radius. 
For the most extended galaxies, the BCG radius and mass 
scale as oc M. For the compact galaxies in the RS09/5 in- 
terpretation thi s scaling is almost re oc . This latter scal- 
ing agrees with Ivan Dokkum et"aLl l|2010h who studied the 
growth of compact galaxies at fixed number density finding 
Te oc Thus for smaller sizes, we get a stronger size 

evolution in better agreement with observation. Seeing that 
our weighting scheme produces central galaxies with appro- 
priate stellar masses for their host dark halo, we now address 
the issue of the change in inner slope of the dark matter pro- 
file in the presence of a baryonic component, and how this 
varies depending on the relative distribution of stars and 
dark matter. 




Figure 3. Normalised radial distribution of stellar mass in the 
z = BCG coming from different contributing galaxies for model 
A (each colour). The solid line shows the stars coming from 
the largest progenitor, taken to be the proto-BCG (in-situ). The 
dashed-dotted line is the radial distribution of stellar mass con- 
tributed by accreted galaxies. The top and bottom panels show 
the results for two interpretations (RS09/1 and RS09/5). In the 
RS09/5 case, accreted material contributes 50 percent of the mass 
budget already at 30 kpc while this occurs at 60kpc for the 
RS09/1 case. The arrows point out the radii containing 10, 25, 50 
and 75 percent of the BCG stars. 



3 EVOLUTION OF THE DARK MATTER 
SLOPE 

3.1 Methodology 

In relaxed clusters, BCGs are usually coincident with the 
centre of the cluster, defined theoretically as the bottom of 
its gravitational potential well and empirically as the centre 
of the X-ray emitting hot gas. Thus we define the centre of 
the galaxy (BCG) and cluster as the position of the parti- 
cle with th e minimum potentia l identified by the Subfind 
algorithm (jSpringel et al.l [20011 ) .We then compute density 
profiles in spherical shells around this centre, using 44 bins 
spaced logarithmically between O.lkpc and 2500kpc. The in- 
trinsic slopes of the density profiles are computed by numer- 
ical differentiation using a 3-point Lagrangian interpolation 
as m I Navarro et al] l|2010l ). 
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3.2 Results for the original RS09 simulations 



We begin by analysing the simulations as originally pre- 
sented in RS09. In Figure 4, we present the intrinsic log- 
arithmic slope 7 = —d\n{p)/d\nr as a function of radius 
r for the dark matter (red dashed lines) and for the total 
matter (solid black lines). The dashed horizontal line 7 = 1 
marks the asymptotic value that the NFW profile should 
reach as r — >■ 0. Results are given for the main BCG pro- 
genitor at 2 = 3 in the left panels and the final BCG at 
z = in the right panels. The upper panels are for run A 
(contracted) and the lower panels for run B (uncontracted). 
The BCG progenitors at z = 3 have inner slopes of 7 ~ 1.3 
and 7 ~ 1 at 5kpc for the dark matter in the contracted 
and uncontracted models, respectively. In the final systems 
at 2 = these latter inner slopes reach 7 ~ 0.9 and 7 ~ 0.8. 
These are depressions of 0.2 < A7 < 0.4 with respect to 
the initial slope at this radius. However, as already noted 
these initial conditions are inconsistent with observations of 
z = 3 galaxies assuming a ACDM universe. The behaviour 
here merely serves to prove that coUisionless evolution can 
reduce the slope of the inner dark matter profile. We still 
need to test whether this holds for galaxies with more re- 
alistic initial configurations, so we turn to results from our 
weighting scheme. 



3.3 Dark matter slope evolution for other BCG 
stellar mass profiles 

Figure 4 shows the initial and final slopes of the dark mat- 
ter density profiles for each of our reinterpretations of the 
two simulations using the weighting scheme presented above. 
The total mass does not of course change from one interpre- 
tation to another, only the partition between stellar and 
dark matter distributions does. Because we have decreased 
the stellar mass by a factor of ten, the dark matter mass in- 
creases to compensate. The distribution of the dark matter 
varies however, according to the concentration we assume for 
the stellar mass. In the RS09/1 interpretation, there is addi- 
tional dark matter at all radii and 7 becomes more negative 
for the dark matter everywhere. In the RS09/5 interpreta- 
tion however, stars are much more concentrated which re- 
sults in a shallower initial slope near the centre. The slopes of 
the dark matter profiles retrieved from the weighting scheme 
are represented by the solid coloured lines. 

The initial dark matter profiles in each interpretation 
differ from each other in the inner regions where stars con- 
tribute significantly to the mass budget (5 — lOkpc) and 
asymptote to the same profiles at larger radii. At z — 3, the 
slopes in model A range from 7 ~ 1.4 to 1.0 at 5 kpc, and 
7 ~ 1.6 to 1.5 at 10 kpc between RS09/1 and RS09/5. At 
z = the final slopes take values 7 ~ 0.5 to 1.0 at 5 kpc and 
7 ~ 1.2 to 1.0 at 10 kpc. These are changes with respect to 
the initial slopes of 0.4 < A7 < 0.5 at both 5 kpc and 10 kpc. 
Note that in a few cases (RS09/4 and RS09/5 model A) the 
final slopes are already shallower than 7 = 1 below 10 kpc, 
a region of interest for observations using stellar kinematics. 

For our model B, at z = 3 the slopes range from 7 ~ 1.3 
to 1.0 to 7 ~ 1.6 to 1.4 at 5 and 10 kpc respectively. At 
z = we also observe significant changes in the inner-slope 
between 0.3 < A7 < 0.4 and 0.4 < A7 < 0.6 at 5 and 
10 kpc. We caution that interpretation RS09/5 has a steep 
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Figure 4. Slopes of the total matter and dark matter profiles at 
z = 3 and z = for the four interpretations and for the original 
RS09 simulations. The top and bottom panels show results for 
models A and B respectively. The dashed red line is the slope of 
the dark matter profile in the original RS09 simulations (prior 
to rescaling). The horizontal dashed black line 7 = 1 marks the 
asymptotic limit for a NFW cusp. The arrows mark the resolution 
limit 3e = 4.1 kpc. The lower coloured arrows mark the half-light 
radii for the different representations. 



slope profile and that for this particular interpretation the 
slope change is much more significant A7 ~ 0.7. 

Thus, we see that the various mergers between galax- 
ies not only change the slopes of the total matter density 
profiles but also those of the dark matter profiles. We show 
that even when baryons account for just one percent of the 
total mass of a galaxy, in agreement with recent abundance 
matching results, their effect is still significant enough to 
have an impact on the scales probed by observations, mak- 
ing initial cusps shallower. In the new interpretations, the 
changes in the slope are more localised to the inner regions 
(5 — 10 kpc) of the galaxy when compared to results from 
the original RS09 simulations where the slope was shallower 
than 7 = 1 out to ~ 30 kpc for model B. 

Aspects of this evolution are pr esent in the explana- 
tion proposed by El-Zant et al.l (|2004h for the shallow slopes 
measured by ISand et aL l|2004h : dynamical friction as the 
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galaxies orbit in the cluster heats up the central dark mat- 
ter and makes cusps shallower. A limitation of this work was 
the use of massive particles to represent galaxies which can 
not capture other processes such as mixing and stripping 
of stars and dark matter. Their shallow cusp was produced 
by evacuating the central dark matter and replacing it with 
massive point-like galaxies. The relative distribution of stars 
and dark matter was thus not followed realistically. 

In our numerical experiment, the galaxies are repre- 
sented by many particles giving them a gravitationally self- 
consistent phase-space distribution. The numbers of mergers 
and of galaxies falling to the centre of the cluster are given 
by the simulation. Using the fact that the galaxies are set up 
initially to be spherically symmetric and in a steady state, it 
follows from Jeans' theorem and Eddington's formula that 
we can construct different equilibrium stellar density distri- 
butions out of the total mass distribution of particles. Our 
simulations show that differences in the weakening of cusps 
from one interpretation to another are due to differences 
in the mixing between stars and dark matter. Within the 
ACDM context, the many coUisionless mergers experienced 
by the BCG produce a weakening of dark matter cusps. 

This is a qualitative result and certain aspects of the 
matter distributions in our simulations remain unrealistic. 
As a result we do not attempt to reproduce observations of 
Abell 383 (e.g. the velocity dispersion and luminosity profiles 
in Nil). Even with a fully realistic simulation such a goal 
would likely be unattainable, given that we have only one 
realisation of a cluster at our disposition. BCG formation 
and evolution is tightly coupled to that of the cluster as a 
whole, so BCGs have a variety of profiles. 

In reality the situation is more complicated, galaxies 
do not form only through coUisionless mergers. Disspational 
processes certainly play an important role in shaping galax- 
ies at earlier times. Baryons dominate the visible regions of 
the progenitor galaxies and can alter the initial distribution 
of dark matter in n umerous ways: through gas expulsion 
I Navarro et "all 1 19961). super nova driven impul sive heating 
i Pontzen fc Governatd l201ll ). AG N feedback llDuffy et al'l 
2010l'l or adiabatic con traction (|Blumenthal et al.l Il986t 
Gnedin et alll2004 12OIII ). 



4 CONCLUSIONS 

We have studied the formation of a 10^^ M© cluster in a 
ACDM Universe, following the evolution of galaxies through 
coUisionless mergers from z — 3 down to 2: = 0. We showed 
that as a result of mixing between stars and dark matter 
in dissipationless mergers, initial cusps (0.8 < 7 < 1.3) can 
be substantially weakened (0.3 < 7 < 0.9) at the inner- 
most resolved radii. Our results indicate that observations 
of shallow dark matter cusps at the centre of clusters are 
not necessarily inconsistent with CDM. We find changes of 
dark matter profile slope at a fixed radius of the order 0.3 < 
A7 < 0.5. 

Another interesting result from this study is that the 
evolution in size is stronger than in mass for more com- 
pact stellar distributions. We find that the trend moves from 
Te <x. M for extended galaxies to Ve oc for more com- 
pact ones. The latter evolutionary tren d has in fact been 
found in a recent observational study by Ivan Dokkum et al.l 



(I2OIO'). Our numerical experiments have some serious limi- 
tations that we hope to improve on in future work. It would 
be interesting to use initial conditions based on the observed 
z = 2 mass-size relation for galaxies combined with stellar to 
dark matter ratios consistent with abundance matching ar- 
guments in order to study the build-up of BCGs and the fate 
of the most compact ellipticals in ACDM more reliably. This 
would allow us to check semi-analytic desc riptions of the 
proces s in the ACDM context such as that of IShankar et al.l 
(|201ll ). We are currently investigating thi s using simulations 
of the Phoenix project (|Gao et al.ll2012l '). 
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